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Abstract 

Techniques for the field-theoretic calculation of a form factor are described 
and applied to a dressed-fermion state of a (3-|-l)-dimensional model Hamil- 
tonian. Discrete light-cone quantization plays the crucial role as the means 
by which Fock-state wave functions are computed. An ultraviolet infinity is 
controlled by Pauli-Villars regularization. 

I. INTRODUCTION 

There has been assembled a sequence of technologies by which one might eventually com- 
pute hadronic form factors directly from quantum chromodynamics (QCD). These include 
the early work by Drell and Yan and by Brodsky and Drell on the relation of form 
factors to Fock-state wave functions in light-cone quantization. The wave functions can be 
calculated, in principle, by the method of discrete light-cone quantization (DLCQ) proposed 
by Pauli and Brodsky P]. Refinements of DLCQ that permit substantive calculations in 
non-super-renormalizable, (3-1-1 )-dimensional field theories have now been tested by Brod- 
sky, Hiller, and McCartor ; a key role is played by Pauli-Villars (P V) regularization , 
implemented through the introduction of PV bosons to the DLCQ Fock-state basis. In the 
following sections, a description is given of how this sequence comes together, and the steps 
are applied to a form factor calculation in the model of Ref . . 



II. FORM FACTORS AND WAVE FUNCTIONS 

Formal expressions for the Dirac and Pauli form factors Fi and F2 of a composite fermion 
can be obtained [|l],3 by considering the current matrix elements 
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(P + Q, a| J+(0)/P+|P, a) = 2Pi(Q2) , (2.1) 
(P + Q, (t| J+(0)/P+|P, -a) = -4a{Q, - zQ,)P2(Q')/2M , 

where = Jq + is the plus component of the current, P+ = + P^ is the hght-cone 
longitudinal momentum of the initial state of mass M, Q is the photon momentum, and a 
is the spin projection along the z axis. By working in light-cone coordinates and in 
the Drell-Yan frame [|l|, with Q = {0,2Q ■ P/P~^,Q±), the form factors can be expressed 
directly in terms of Fock-sector wave functions ip^"^^ as 

Fi{Q') = E„ Ej e, / 5(1 - E n 

i i 

F2{Q^) = En e,/ 5(1 - E n (2-2) 

i i 



2M 



Here n is the number of constituents; Cj is the charge of the struck constituent; the Xi 
are longitudinal momentum fractions pf /P~^ for constituents with momenta pi] the k±i are 
relative transverse momenta; and 

are transverse momenta relative to the new P + Q direction. Thus, given the wave functions 
one can compute the form factors. 
In light-cone quantization these wave functions can be found by diagonalizing the mass- 
squared operator P+P~ — P]_ = i^LC, traditionally called the light-cone Hamiltonian [§. 
This is simplest in a frame where the total transverse momentum is zero and in a basis where 
P"*" is diagonal. Then only the light-cone energy operator P~ has a nontrivial representation, 
and the eigenstate can be expanded explicitly in terms of momentum Fock states \n : Xi , k±i) 
and the desired wave functions: 

1^, ^) = E / 5(1 - E ^0 n ^16vr3'^(E k^^) n 1^ ^ ^^^) • (2-4) 



It is a significant advantage of light-cone coordinates that such an expansion is well-defined, 
in the sense that there are no contributions from disconnected vacuum pieces P]. Another 
advantage is that the boost invariance of x and k± permit the same wave functions to be 
used in the construction of the boosted state \P + Q , a) . 

III. COMPUTATION OF WAVE FUNCTIONS 

The wave functions can be computed by applying DLCQ to the eigenvalue prob- 
lem ifLc|-P)Cr) = M^\P,a). Pauli and Brodsky applied this method to (l-l-l)-dimensional 
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models for its initial trials 0. Since then much work has been done [§, particularly in 1+1 
dimensions. The greater complexity of (3+l)-dimensional theories has made progress there 
much slower. The need for regularization and nonperturbative renormalization has been 
most telling 0. 

Recently a scheme for use of Pauli-Villars (PV) ultraviolet regularization within 
DLCQ calculations has been successfully tested for (3+l)-dimensional model Hamiltonians 
IQJ^. Massive PV bosons are introduced as additional constituents in the Fock basis, with 
imaginary couplings chosen to produce desired cancellations in perturbation theory. The 
eigenvalue problem is solved nonperturbatively, and the limit of infinite PV mass is taken. 

Of the two model Hamiltonians considered, the more sophisticated is |^ 



+ 



167r^g+ 
dpi(Pp± 



+ gi t , fj'l + ql t 

q+ / P+ "^i^i q+/P+ "'H^'^l 
dp2d'^pj_2 f dq^d'^qx^ 



f dpjd^p±i r dj9^d"p_L2 f dq-^d^q^ ^ ^ 



lQ'n-^P2 



(3.1) 



X 



+ia\g6{p^ -P2 + q) + iaiqKli ~ 2.2 



q) 



and spin a, act as 
The fermion mass 



Fermions created by 6],^., with light-cone momentum p = {p'^,Px,Py] 

sources and sinks for physical and PV bosons created by and a|g 
is M, and the boson masses /i and /ii. The imaginary coupling of the PV boson causes a 
cancellation of an infinity in the fermion self-energy. The M^p^/P^ counterterm is then all 
that is needed to remove the shift in the mass, with Mq ~ In fii/fi. 

The eigenvalue problem for this Hamiltonian, in the one-fermion sector with total mo- 
mentum P, reduces to a system of integral equations 



P+/P+ 



my IP 

dq^d'^q± 



E 



E 



III 



rj/P+ 



v/lGvr-^g 



(n+l,ni) 



{q,,q,rj,p- q) 



n 



(n— l,ni) 



+ i + 1 y 



lOvT'^g/ 
dr^di^r^ 



,q„,rj,p + q. 



+ 



1 



(7i,ni+l) 



{q^,Lj,L,p-id 
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(n,ni — 1) 



'-i-i'-i+i' 



-n\ 1 



p + r •) 



(3.2) 



for the wave functions where n and n\ are the numbers of physical and PV bosons. 

There are two bare parameters Mq and which are fixed by setting values for physical 
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quantities chosen to be M and (:0^(O):) = (-P, cr| :0^(O): |-P, cr). The latter was chosen for 
ease of computation, in the form 



= E ny dqtd'q^^ll J dr^dh^, (3.3) 

1 n j 

2 



n=l,rai=0 i 




These renormahzation conditions must be solved simultaneously with the eigenvalue prob- 
lem. 

The DLCQ method translates the field-theoretic eigenvalue problem into a matrix prob- 
lem. Discrete momentum values p"*" = / L and p± = {nxTr/L±,ny'7i/ L±) are used, with 
L and L± the length scales associated with the approximation. Momentum integrals are 
approximated by trapezoidal sums over the discrete points. Near the endpoints special 
weighting factors are needed for better accuracy, as discussed in Ref. 

The length scales L and L_i_ define a coordinate-space box within which the fields are 
assigned periodic boundary conditions, except for a longitudinal antiperiodic condition for 
the fermion. The integer n is then odd for fermions but even for bosons. The total longi- 
tudinal momentum for a one-fermion state defines an (odd) integer K = P+L/vr called the 
harmonic resolution |^. The transverse integers and ny range between — A^_|_ and N±, 
with A''^ fixed by a cutoff + p]_- < A'^pf /P~^. The cutoff is needed to produce a matrix 
of finite size but is not used as a regulator. The parameters K, L±, and determine the 
numerical approximation. 

The matrix eigenvalue problem is readily solved for the lowest massive state by use of the 
complex-symmetric Lanczos algorithm [1^. The algorithm allows the matrix to be stored 
in a compact form and to be referenced only in matrix-vector multiplications. This feature, 
combined with the extreme sparseness of the matrix, has permitted calculations with as 
many as 10.5 million basis states f^. 



IV. COMPUTATION OF A FORM FACTOR 

With these techniques in place we may consider calculation of Fi for the dressed fermion 
state of the model.Q In Ref. only the slope -Fi'(O) was computed, from an expression 
derived from ( |2.2|) . Here we compute with ( p.2|) directly. 

The structure of the Fi formula is that of an overlap of momentum wave functions, with 
one shifted in a transverse direction, which we take to be x so that = Q±x. The shape 
of the one-boson wave function is shown in Fig. |l|. For the coupling strength considered, 
this wave function provides the primary contribution to the variable part of Fi. Of course, 
the bare fermion Fock state provides a Q_L-independent contribution equal to that state's 
probability. 



The interactions of the model do not flip the spin, and F2 is therefore identically zero. 
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FIG. 1. A cross section of the boson- fermion two-body amplitude taken at fixed longitudinal 
momentum fraction x = 5/9 and at fixed ky = 0, with K = 9, {:(j)'^{0):) = 1, and ^1 = 
The cutoff A-^ and the transverse resolution N± are varied to keep the transverse scale L± fixed 
at one of the following values: 1- (black), \/2- (gray), and 2- (white). Different symbol shapes 
correspond to different values of A^j_. The peaks are normalized to be equal at kx = 0. The points 
at zero amplitude mark the transverse range, which is set by the cutoff. 



Calculations have shown that the wave functions quickly become independent of K 
and L± as these parameters are increased. The independence with respect to L± can be seen 
in Fig. |l|. However, -F'(O) was found to be sensitive to L^, and one would expect the tail of 
Fi to be sensitive to A^ as well as L±. That this is the case can be seen in Figs. ^ and ^. For 
small A^ the form factor quickly reaches the bare-fermion contribution. For larger A^ the 
bare-fermion contribution increases slightly, and the approach of Fi to this limiting value 
becomes more gradual. When A^ is large enough for the shape of Fi to appear converged 
there remains a significant L± dependence, as seen in Fig. |^, even though the one-boson 
wave function changes little with L^. The form factor remains sensitive because L± controls 
the approximation to the integral of the wave function product from which Fi is computed. 



V. SUMMARY 

These results show that, at least for model Hamiltonians, a field-theoretic calculation 
of Fock-state wave functions and bound-state form factors can be carried out. The added 
Pauli-Villars particles provide the ultraviolet regularization without making the basis size 
unmanageable. Work on a more complete field theory, a single-fermion truncation of Yukawa 
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FIG. 2. The form factor Fi for fixed resolution K = 9 and L_i_ = 2t: / ^. Various cutoffs are 
considered. The model parameter values are (:(/>^(0):) = 1, = jj?, and ii\ = 

theory, is underway. Consideration of the dressed-electron and positronium states of quan- 
tum electrodynamics would be a natural next step. QCD will require a more sophisticated 
approach, perhaps relying on heavy supersymmetric partners to play the role of the Pauli- 
Villars particles. 
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FIG. 3. The form factor Fi for fixed longitudinal resolution K = 9 and varying transverse 
resolution L± . The largest available cutoff was used in each case. The model parameter values are 
(:</,2(0):) = 1, m2 = /x^, and /xf = lO/i^. 
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